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Abstract 

We present a study of the recently discovered spatially-extended chaotic 
state known as spiral-defect chaos, which occurs in low-Prandtl-number, large- 
aspect-ratio Rayleigh-Benard convection. We employ the modulus squared of 
the space-time Fourier transform of time series of two-dimensional shadow- 
graph images to construct the structure factor S{k,u;). This analysis is used 
to characterize the average spatial and temporal scales of the chaotic state. 
We find that the correlation length and time can be described by power-law 
dependences on the reduced Rayleigh number e. These power laws have as 
yet no theoretical explanation. 
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I. INTRODUCTION 



Attempts to understand the complex deterministic motion of spatially extended nonlinear 
dissipative systems, so-called spatio-temporal chaos (STC), are at the forefront of research in 
nonlinear dynamics . In contrast to simple chaotic systems in the time domain, in which 
a few important degrees of freedom are nonlinearly coupled, spatially-extended systems 
seem to require an infinite number of modes for their description. Simple chaotic systems 
are typically described by a few amplitudes which are determined by coupled nonlinear 
ordinary differential equations. Extended systems, on the other hand, must be described by 
nonlinear partial differential equations for which, in general, no natural reduction to a few 
modes is known. This extended form of chaos has been characterized to some extent for 
a variety of experiments and simulations, but continues to resist any understanding based 
on generalizations of the few-modes case. In this paper, we describe experiments on the 
spiral-defect-chaos (SDC) state in Rayleigh-Benard convection |,|,|ig (RBC). RBC 
occurs beyond the first hydrodynamic instability of a thin horizontal fluid layer heated from 
below, and is particularly suitable for the study of complex pattern-formation phenomena 
because so much of the theoretical foundation has been established by the extensive stability 



analyses of Busse and coworkers. |Tl|-p!^ 

The SDC state consists of a disordered pattern of convection rolls, which prominently 
features the persistent spontaneous appearance and disappearance of rotating spiral defects. 
It offers an unusually controllable experimental system for the study of STC. Other examples 
of STC which have been studied experimentally include for instance "dispersive chaos" in 



binary fluid convection |]T5|, chaotic regimes in electro-convection in nematic liquid crystals 
HTgHTsI , chemical-reaction patterns , parametrically-excited surface waves pO| , and RBC 
under rotation in the Kiippers-Lortz unstable regime Most of these cases have the 

common feature that STC is encountered following the breakdown of an ordered pattern as 
a control parameter is increased. The onset of STC in such systems is typically linked to 
the appearance and proliferation of defects in an underlying pattern, leading to complex, 
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persistent pattern dynamics. These patterns, although disordered, still have a characteristic 
wavenumber and are far simpler than fully developed turbulence in which a broad spectrum 
of length scales is involved. Usually, a crucial role is played by the size of the system, which 
must be large when measured in units of the characteristic size of the pattern elements. 
Only in the limit of large system size, may one hope to approach the limit of a statistically 
homogeneous chaotic state that is free of boundary effects. 

RBC may be parameterized by four dimensionless numbers. The Rayleigh number R is 
a dimensionless temperature difference across the layer and is given by 

B.^?^^. (1) 

It has a critical value of Rc = 1708 for a laterally infinite layer at the onset of convection. 
Here AT is the temperature difference, d the layer thickness, g the acceleration due to gravity, 
a the thermal expansion coefficient, k the thermal diffusivity and u the kinematic viscosity. 
It is convenient to define the reduced Rayleigh number e = {R/ Rc) — 1, so that e = 
corresponds to the onset of convection. This number may be regarded as the experimental 
control parameter. The Prandtl number, given by cr = z//k, is a second parameter that 
describes the fluid properties. It has important effects on what secondary instabilities are 



expected in the pattern that is found above onset |ll3|Jl^ . In our experiments, the working 
fluid was pressurized CO2 gas, for which o" ~ 1. This value of a is characteristic of gases. 
A Prandtl number close to or less than one seems crucial for the appearance of SDC The 
lateral extent of the convecting layer is described by the aspect ratio T = r/d, where r is 
taken to be the radius for a circular geometry. The SDC state is found only for F ^ 1. 
The lateral boundaries may also influence the pattern via roll-orienting effects which are not 
completely specified by F alone. We shall show that these effects are unimportant to SDC, 
provided F is sufficiently large. Finally, the non-Oberbeck-Boussinsesq (non-OB) parameter 
V, defined in the next section, characterizes the degree to which the fluid properties 
vary over the height of the layer. Provided this parameter is small enough, as will be the 



case for the experimental conditions we consider here, the Boussinesq approximation ||T0 



applies, and the pattern near onset consists essentially of rolls. 

SDC was first observed in gas convection using CO2 under non-OB conditions, that 



is when V was not small. 22 In this situation, the onset of convection is a transcritical 



bifurcation to a pattern of hexagonal flow cells p2| , p3| . This pattern became unstable at 
e ~ 0.1 to a roll pattern which, at higher e, became the disordered time dependent SDC 
state. In later experiments it was found that SDC persisted under OB conditions. 
These conditions can be achieved experimentally at higher gas pressures 0, or in thicker 
samples P, ^ , p5| . In the experiments with circular cells discussed here, for which F = 74.6, 



SDC begins to appear near e ^ 0.25 and is well developed by e ~ 0.5. We found that this 
state is robust in the sense that it is found in both circular and square cells, independent of 
the temporal history of e, and with various lateral boundary conditions. Similar experiments 
have also identified this state in circular cells with both smaller and larger PJ26[ 



aspect ratios. SDC has been observed using several other pure gases with Prandtl numbers 
near one, including SFg [p|,p5[, N2 [^], and Ar [^. In experiments near the critical point of 
SFg, Assenheimer and Steinberg found that SDC evolves into a state of "target chaos" 
as the Prandtl number is increased. SDC also has been studied in experiments in which 
there is an additional control parameter, namely the rotation rate of the convection cell 
about a vertical axis [§. It also has been observed recently in various gas mixtures ||25 



which have smaller Prandtl numbers than pure gases. SDC has been found in several 
numerical simulations of convection patterns, as modelled by modified Swift-Hohenberg 



equations p9|-|3T|. Decker et al. ||3^ were able to quantitatively reproduce the SDC state 
observed in our experiment by a Galerkin truncation of the solutions to the Navier-Stokes 
equations. 

The apparently ubiquitous nature of SDC may seem at first to be in conflict with the 
well-known picture developed by Busse and Clever |]T3| , p!^ of the secondary instabilities 
available to straight convection rolls. However, their analysis tested the stability of straight 
roll solutions of the convection equations to arbitrary perturbations. The SDC state, which 
involves defects and strong roll curvature in an essential way, appears to correspond to a 



different attractor that coexists with that of straight rolls over a range of parameters. It 
is apparent from our experiments with square cells, discussed below, and from previous 
studies by Croquette ||33|, that defect-free straight rolls represent a rather special situation 



that is extremely difficult to prepare in large-F cells without special forcing schemes. With 
small-amplitude random initial conditions, the SDC attractor is always found in simulations 



3J|. In the experiment, random fluctuations which nucleate the growth of convection 



patterns, or the roll curvature induced by sidewalls in round cells or by defects and grain 
boundaries in square cells, invariably seem to put the system into the attractor basin of SDC. 
We remark that one of the outstanding problems in the study of STC is how to generalize 
the concept of "basin of attraction", which here can only be understood rather loosely in 
analogy with low-dimensional chaos. A direct comparison of the SDC state with the stability 
boundaries given by the Busse-Clever analysis has been presented previously PJ^l- 

In the absence of a general framework for the understanding of STC, we fall back on 
statistical descriptions of the SDC state. There is however, every reason to believe that SDC 
is deterministic in the sense that the external noise driving the system is quite small 



and that the observed fluctuations represent the deterministic outcome of complex dynamics 
intrinsic to the system. In order to characterize the statistical structure of SDC, we made 
use of the space-time structure of the state, in the form of time series of snapshots of the 
pattern. We analyzed these using the three-dimensional structure factor S{k,uj), which we 
deflne to be the modulus squared of the Fourier transform of the space-time shadowgraph 
data. Many other sorts of analysis are possible, but this method has the virtue that the 
results may be interpreted in a nearly model-independent way to extract the global spatial 
and temporal scales of the dynamics. This is analogous to familiar and standard techniques 
for the study of disorder in materials by, for example, scattering experiments. A great deal 
of complexity is suppressed by this procedure, which amounts to a massive averaging over 
the chaotic fluctuations. The central and open problem of STC remains to understand this 
complexity by a more speciflc reduced description. 

The rest of this paper is organized as follows; in Sec. || we describe the apparatus. In 



Sec. |TT| we describe in several subsections the phenomenology of spiral defect chaos, first 
qualitatively, and then quantitatively using the structure factor. In these sections, we also 
discuss the main results on the spatial and temporal scales of the chaotic state. Sec. fTV] 
presents a discussion of SDC in the context of the general features expected of STC Finally, 
Sec. |V] is a short conclusion. 



II. THE EXPERIMENTAL APPARATUS 

The experimental apparatus was one of several of similar design that have been used in 
other experiments |]3|,^,|6|-p|,^IH26|l . A comprehensive description has been given here we 



describe only some essential features. The top surface of the convection cell was a circular 
9.5 mm thick, optically fiat sapphire window. The window was temperature controlled 
by circulating water held at the same pressure as the convecting gas. This avoided the 
deformation of the window that would result from a pressure differential across it, and made 
cells with an extremely uniform thickness possible. The residual distortion in the flatness of 
the sapphire, which was ~ l/im, was due to stress from its mechanical mounting. A pressure 
vessel enclosed the cell and a closed volume of water which was circulated over the top plate 
by a small pump. The temperature of the circulating water was fixed at 24.00 ±0.02°C and 
regulated by a computer controlled servo. The water pumped across the top sapphire could 
be regulated to ±0.2mK. 

The bottom surface of the cell was an aluminum plate, 9.5 mm thick, which was dia- 
mond machined to ~ 1/im surface flatness. The plate was heated by a commercial thin film 
heater which covered its entire bottom surface. Its temperature was sensed by an embedded 
thermistor. The temperature of the plate, which was under the control of a computer servo 
similar to that used for the circulating water, was the main experimental control parame- 
ter. It varied between about 30°C and 45°C, measured to ±0.02°C, and was regulated to 
±0.2mK. The position of the bottom plate was determined by an adjustable support which 
included three piezoelectric legs. These allowed adjustment of the cell thickness over a range 
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of ~ 10/im, and were used as a fine adjustment after the cell was pressurized. The height 
of the cell and its uniformity were measured interferometrically at the working pressure, as 
described in Ref. I^Gj. The cell-height uniformity was typically ±l/im. The cell thickness 
could be determined to within ±lyum at a reference bottom plate temperature. As the 
bottom plate temperature was increased, the cell thickness changed, due primarily to the 
thermal expansion of the aluminum bottom plate. This expansion was about 1% of d at the 
highest bottom-plate temperatures. It was measured, and a correction was applied in the 
calculation of e described below. 

The lateral sidewalls of the cell were defined by a porous gasket made of three sheets 
of filter paper. The relatively low thermal conductivity of the paper reduced temperature 
gradients caused by the tendency of the gasket to act as a thermal short circuit between the 
top and bottom plates. The effect of such gradients on patterns in circular cells at low e is to 
stabilize cell filling spirals p3| , p7| or concentric target patterns []38| , |39| , |33| ,^,|37| in which the 



rolls lie parallel to the sidewalk We constructed both circular and square cells using paper 
sidewalls for which the gradients were sufficiently small that the rolls lay perpendicular to 
the sidewalls except at very small e. In order to study the effect of sidewall forcing, we 
embedded fine heater wires in the paper sidewall which were sufficient to restore strong 
parallel roll alignment near onset. In the circular cell, the sidewall heater extended around 
the entire circumference, while in the square cell there were heaters on two opposite edges. 
The circular cell was 43.88 ± 0.05 mm in radius and 588 ± 1/im high, giving an aspect ratio 
r = radius/height ~ 74.6. The square cell was 62.5 ±0.1 mm on each side and 639 ± 2/im 
high, so that the aspect ratio V = width/height = 98. In both the round and the square cells, 
the CO2 gas was held at a constant pressure of 33.1 ± 0.1 bar. The pressure was regulated 
to be constant to within ±0.01% by servo controlling a heater in a gas filled ballast volume 
connected to the cell. 

The convection patterns were imaged using the shadowgraph technique [^,0 which 
maps the deviation of a parallel beam of light reflected off the mirrored bottom plate. The 
deviations are caused by the varying index of refraction of the gas due to the small lateral 



temperature differences across a convection roll. In these images, the white areas correspond 
to cool downward ffowing gas, while the dark areas are warm upffows. To make quantitative 
pattern measurements, the shadowgraph was operated in a linear regime and the image was 
captured by a CCD video camera using a computer and an 8 bit frame grabber. Dividing 
the image by a reference image taken below the onset of convection largely compensated 
for nonuniformities in the illumination. Our quantitative results are based directly on these 
ratio images. For the purpose of visual presentation, the ratio images were re-scaled so as 
to render them with increased contrast. 

Values of the relevant dimensionless parameters were calculated from the known gas 



properties. |Q The properties were evaluated at the average temperature T = {Tf + Tb)/2 
where Tt and Tf, are the temperature at the top and bottom of the fluid layer respectively. 
The temperature was that measured by the thermistor embedded in the bottom plate. 
The temperature Tt was slightly larger than that of the circulating water because of the 
finite thermal conductivity of the sapphire top plate, which was taken into account. At the 
onset of convection this difference was 0.06°C. We calculated the vertical thermal diffusion 
time scale t„ = d"^/ /t as a function of T and used it to rescale all the times in the experiment. 



The non-OB parameter V is given by |T2 
P = AT 



Pi f d^p\ P2 /dri\ P3fdA\ p^ (dCp 



T=T 



(2) 



where a is the isobaric thermal expansion coefficient —{1/ p){dp/dT)p, A the thermal con- 
ductivity, and Cp the heat capacity at constant pressure. The coefficients pi are dimen- 
sionless linear functions of and are given in Ref. The derivatives in Eq. (0) were 



evaluated numerically from the fits to p, rj, A and Cp given in Ref. Since only Tb 



was changed while the temperature of the circulating water bath was held fixed, the mean 
temperature T changed and consequently a, and V all varied slightly with e. 

In the circular cell, we found the onset of convection at Ton = 27.32°C and 
AT/^P* = 6.513±0.005°C. The result for AT/^^* agrees to within about 3% with the value 
ATJ^^ (Ton) = 6.74°C calculated from the gas properties and the thickness. In order to 
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calculate e well above onset, we must take into account the variation of the fluid parame- 
ters, and hence AT^, with T. We define e = [AT/ATco(T)] — 1, where ATco(T) is the critical 
temperature difference for a fluid with properties appropriate to the mean temperature T. 
To remove the effect of the small systematic difference between ATc^^^ and AT J^* {Ton), we 
defined ATco(T) as 

Areo(f ) = AT/**(f ) + [AT/^P* - ATj'\fon)]. (3) 

Thus, we use fits to the gas properties to get the functional form of ATc-^**(f ), which is 
nearly linear, and include the small correction term in square brackets in Eq. (j^) so that this 
curve passes through the experimentally observed value of ATc*^^^*. This correction becomes 
relatively unimportant as e becomes large. Using this definition of e, the control parameter 
dependence of V in the circular cell is given approximately by P ^ —0.94 — 1.22e, while the 
mean Prandtl number was cr = 0.94, and the vertical thermal diffusion time was = 1.44s. 
For our range of AT, a varied by 3% and by 10% due to the variation of T. 



III. SPIRAL DEFECT CHAOS 



A. Qualitative features 



In this section, we describe the general features of the SDC state, and how they depend 
on some simple variations of the experimental conditions. MBCA96Figures |l| and |^ shows 
the sequence of patterns observed in a circular cell with the sidewall heater turned off. For 
small e ^ 0.05, as shown by Fig. |l|a, the pattern consists of motionless straight rolls, with 
small grain boundaries around the circumference. This is consistent with the prediction 
of Schliiter, Lortz, and Busse |Tl[] and the weak forcing nature of the paper sidewalls. As 



e is increased, the tendency of the sidewalls to enforce a perpendicular roll orientation 
becomes more pronounced, eventually leading to the frustrated pattern shown by Fig. |l]b 
at e ~ 0.245. These patterns typically show at least three focus singularities around their 
circumference and several slow-moving defects and grain boundaries in the interior. The foci 



often act as sources or sinks for rolls and dislocations This state may be contrasted 
with the two-focus "Pan-Am" pattern seen in smaller F cells p9| . |33| J^ in this regime. At 
e ~ 0.25, the first spiral defects are seen. These appear quite intermittently, and usually 
develop within highly defected regions near the center of the cell. The spirals exhibit a 
tendency to rotate and evolve on shorter time scales than the rest of the pattern. By 
e ~ 0.5, the nucleation of spiral defects is continuous, usually in the interior of the cell, 
and the highly defected SDC regions coexist with slow-moving focus singularities around 
the edge of the cell, as shown by the Fig. ^ja. As e is increased, the SDC state progressively 
invades the entire cell, as shown by Fig. ^3. Even at the highest e ~ 1.4, the rolls continue 
to meet the sidewall perpendicularly, creating a circumferential fringe of small focus defects 
as long as the sidewall heater is turned off. 

We may note several qualitative features of the SDC state, as it appears in the interior 
of the cell, away from boundary effects. First, we find that both right and left handed 
spirals appear, in roughly equal numbers, and are not created in pairs. In fact, the spiral 
cores are only the most visually striking defect in a pattern which contains many other 
defects including dislocations, grain boundaries etc. One-armed spirals are certainly the 
most common, but target patterns and two-armed spirals also appear occasionally. The 
processes by which spirals are created and destroyed are not simple in general, but some 
common scenarios have been described by Assenheimer and Steinberg The rotation 
rate of the spirals is variable and difficult to estimate except in a few particularly long- 
lived examples. In cases where a sustained rotation is seen, the spiral core rotates in a 
"winding-up" sense, i. e. so that with successive turns rolls move outward from the core. 

The bottom plate temperature could either be ramped slowly, or stepped rapidly. We 
are only concerned here with patterns which have reached a statistical steady state. In the 
absence of other information, a pattern transient may be expected to take at least several 
horizontal thermal diffusion times th = F^t^, ^ 2.2h. In fact, over most of the range of e we 
studied, the SDC state takes much less time than this to establish itself, requiring only a 
few hundred t^. This is apparently due to the surprisingly fast time scale of SDC, as we 
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will make quantitative below. Nevertheless, our experimental protocol was to step e rapidly 
(within ^ lOt^) to a target value, then to wait at least 2th before collecting data. Runs in 
which e was ramped slowly, (ti,^ — 10^^), produced similar patterns at low e and a similar 
transition to SDC 

In order to examine the effect of changing the sidewall boundary conditions, we performed 
a few runs with the sidewall heater on. For small e ^ 0.20, the main effect is to convert 
straight roll and frustrated patterns similar to those of Fig. ^ to large cell-filling spirals, 
as shown in Fig. ^. These large spirals were qualitatively similar to the rotating spirals 



observed previously P3| , p7| . If the sidewall heater is turned off after a cell- filling spiral is 
established, it was found that the center of the spiral soon moves to one side and a frustrated 
pattern similar to that shown in the Fig. |l]b develops after a long transient. At higher e, 
with the sidewall heater on, a similar SDC state is found, as shown in Fig. ^. Here SDC 
evolves in the interior of the cell, while at the perimeter the rolls lie parallel to the sidewalls 
and no small foci are present. We have not yet undertaken a detailed study of the transition 



between the large spiral states p3| , p7| and SDC. It would be interesting to compare the 
onset of SDC in a circular cell with and without sidewall forcing, and in various aspect 
ratios, using the methods of Ref. . 



To test the sensitivity of the SDC to the sidewall geometry, we examined the phenomenon 
in our square cell. Figure || shows the patterns observed. We attempted to use the sidewall 
heaters on opposite edges to prepare a state of straight rolls, but this proved to be impossible 
even for very slow ramping of e through onset. This was apparently due to a very narrow 
range of hexagonal flow cells that appears just above onset, due to non-OB effects; these 
tend to decay to roll patches that meet at 120°, leading to persistent grain boundaries as in 
Fig ^(a). We also tried deliberately breaking the symmetry by putting a large thickness (and 
hence e) wedge across the cell from edge to edge, and by tilting the entire apparatus several 
degrees toward one edge. None of these methods were sufficient to prepare straight rolls. 
While this effort would probably be successful with a sufficiently drastic forcing scheme, it 
seems fair to say that defect-free straight rolls are practically inaccessible given the initial 
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conditions and boundary conditions that were experimentally available in our cell. They have 
been successfully prepared in smaller aspect ratio rectangular cells using sidewall heaters, 
however |Q. At higher e ~ 0.5, we found the familiar SDC state in the square cell, 
regardless of the sidewall heating, as shown in Fig. |^(b). 

The observation of SDC with and without sidewall forcing in circular cells, in cells of 
various F, 0] in sufficiently large square cells, and in Navier-Stokes simulations with periodic 
boundary conditions lends support to the conclusion that SDC is a generic state for low- 
er and large-F convection, which is not related to the detailed boundary conditions. In the 
next section, we take up the quantitative characterization of SDC. This will show that the 
characteristic length scales in the SDC regime are much smaller than the size of the cell, as 
one expects for an "extensive" chaotic state. 



B. Quantitative analysis using the structure factor. 

In this section we describe the statistical methodology we used to quantitatively charac- 
terize the SDC. A time series of two-dimensional shadowgraph snapshots is taken and may 
be regarded as a sample of three dimensional space-time. We presume that the states for 
which we take such series are statistically stationary. Our basic tool is the structure factor 
S{k,uj) which we determine experimentally by taking the modulus squared of the discrete 
Fourier transform of the three-dimensional space-time data set. Since phase information is 
destroyed by taking the modulus, this already represents an enormous simplification of the 
data. We further simplify the structure factor using certain windowing and averaging tech- 
niques described below, to extract a few global features of the motion, namely correlation 
lengths and times, as a function of e. This analysis leaves much to be desired and averages 
out a lot of the interesting details of the dynamics; it has the virtue of simplicity and is 
almost completely free of ad hoc assumptions. 

Two types of spacetime data were collected. In the first type, time series were taken 
with a long wait, typically several hundred ty, between images. This time step was chosen 
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to be long compared to the typical timescale for pattern change, so that successive images 
were nearly uncorrelated. This type of data was used for obtaining long time averages of 
S{k, t). A second type of time series consisted of images taken a few apart, so as to obtain 
temporal correlation information. In both cases, each divided image was 256 x 256 pixels 
with 8-bit greyscale resolution. One spacetime sample usually contained 256 such images. 
The size of the square image was such that its corners just met the edges of the round cell. 
In order to avoid including the pattern in the corners of the image, which is influenced by 
boundary effects, and to reduce the aliasing effect of the sharp edged picture, we prefiltered 
the images by multiplying them by a radial Hanning function given by 

H{r) = [1 + cos{nr/ro)]/2, r < Tq (4) 

H{r) = 0, r> ro. (5) 

The radius ro was O.TIF in units of d. A one-dimensional Hanning window was also 
applied to the spacetime data in the time direction. The effect of each of these windows is 
to convolve the windowing function with the Fourier transform in such a way that ringing 
due to the sharp boundaries of the dataset is suppressed. 

C. The time averaged structure factor 

— * 

We consider the time averaged structure factor S{k) for the temporally uncorrelated 
spacetime data first 0. This could either be obtained by numerically integrating S{k,u) 
over uj, or more simply by time averaging the instantaneous spatial structure factors S{k, t)at 
each time step. The latter is equivalent to the former. When SDC is well developed, S{k) 
shows a broad ring centered on = 0. Only data within an annular band near this ring were 
further analyzed; we ignored, for example, a large spike near k = 0, which is due mostly to 
large scale nonuniformities of illumination. A diffuse positive background, presumably due 
to camera noise, pervades fc-space, and underlies the ring. The shadowgraph was operated 
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in as nearly linear a regime as possible, but at the highest e, a trace of a second harmonic 
ring was detectable. This ring was either due to residual nonlinearity in the shadowgraph, 
or possibly higher order components within each roll of the flow pattern itself. In any case, 
only data near the primary ring were considered, and the very small tail due to the second 
harmonic was treated as part of the diffuse background. 

In the vicinity of the onset of SDC, for e ^ 0.3, the ring in S{k) appears irregular 
or broken, while near the onset of convection it is reduced to two spots. For e ^ 0.4, 



where spirals appear continuously, S{k) was nearly azimuthally symmetric. Hu et al. ||24 
have shown that the uniformity of the power around the ring is a useful diagnostic for the 
appearance of spirals and hence of SDC. This is obviously due to the fact that spirals or 
targets naturally contain rolls of all orientations, while the frustrated patterns similar that in 
Fig. |I]b usually contain a dominant roll orientation. More precisely, nonuniform /c-space rings 
result when the slow roll rearrangements in the frustrated patterns have insufficient time 
to sample all roll orientations within the averaging time of the dataset. Once the quickly- 
evolving spirals begin to reliably appear, the ring fills in over relatively short averaging 
times. 

The inverse width of the ring is a quantitative measure of the length scale of the patches of 
rolls in the pattern. We performed an azimuthal average in /c-space to obtain the azimuthally 
and time averaged structure factor S{k), where k = \k\. Using the average wavenumber 
(k) from this distribution, defined by 

^ j\k\s{k)d'k _ s^es{k)dk 

^ ^ ~ JS{k)d^k j^kS{k)dk ' ^ > 

we define the moments /i„ of S about {k) as 

7(1^1 - {k)YS{k)d^k'\ - {k)YkS{k)dk 



JS{k)d^k 



S^kS{k)dk 



(7) 



The correlation length ^ is calculated from the second moment via 

e = i/v^ • (8) 
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The values of higher moments extracted from data are very sensitive to noise in the wings 
of S{k). Rehable values up to the fourth moment were extracted from S{k) by first fitting 
kS{k) to a smooth phenomenological model function, given by 

kS{k) = [ao + aik] + [ag + a3A;]e-["5(^-''*)'+"«(^-'^*)'+"^('^-'^*)*l , (9) 

using a nonlinear least squares algorithm, where the are adjustable parameters. The ao 
and Oi terms of the model function are a linear background contribution, while the rest 
form a skewed gaussian "bump" piece which decreased exponentially in its wings. With 
eight parameters, an essentially perfect fit could always be achieved. Once the data were fit, 
moments were extracted numerically from the model function, using only the second term 
in Eq. ^ This method was robust to noise in the wings, modifications of the fitting range, 
and omissions of points in the data set. It removed the background contribution as well. 

Figure ^ shows the structure factor that results from this averaging procedure for several 
values of e. The solid lines are the fits of Eq. ^ to the data. For the lowest e shown, the rolls 
are essentially perfectly straight over the whole region sampled. Nevertheless, S{k) still has 
some width due to the finite number of rolls in the sample, and the effect of the Hanning 
window. This may be regarded as a "resolution-limited" structure factor; it corresponds to 
our maximum resolvable correlation length ^ ^ 14d. As e is increased, (k) shifts to smaller k, 
and S{k) becomes significantly skewed to higher k in the SDC regime, as we discuss below. 

In a previous publication [Q, we compared the mean wavevector (k) of the pattern to 
the band of stable wavevectors expected for straight rolls, according to Busse and Clever 
T3| , p^ . The trend in (k) shows no strong changes near the onset of SDC, and roughly tracks 



the middle of the predicted stable band, starting from onset. This may be contrasted with 
the results of Hu et al, who found significant deviations from the middle of the band in a 
r = 40 cell. Evidently, the processes that select {k) at various e depend on F. 

The correlation length ^ decreases with increasing e 0, and our results can be fit by 
the power law ^ oc e-o-43±o.o5_ show that the exponent of this power law is consistent 
with 1/2, we plot vs. e in Fig. 0. If ^ is exactly proportional to e~^/^, this will be a 
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straight line passing through zero. As Fig. ^ shows, this is quite consistent with the data, 
with small systematic deviations only at larger e. A fit of = ^cT^e to the data at small e 
is shown by the solid line in Fig. 0, and gives — 2.3 (in units of d). We have at present 
no theoretical insight into why this power law should appear, why the exponent should be 
so close to 1/2, and why ^ should diverge at e = rather than, say, at some small positive 
e. The correlation length is presumably unrelated to the one which appears in an amplitude 
equation valid near e = 0. That correlation length has the same exponent 1/2, but has 
^0 = 0.385 which is much smaller than our result. 

The qualitative meaning of the correlation length as it relates to the interpretation of 



images requires some clarification It is tempting to assume that ^ is a measure of the 
typical spiral diameter, but this is incorrect. A disordered pattern consisting of a patchwork 
of nearly ordered subpatterns of typical size i will only exhibit a correlation length ^ ^ £ if 
the patches consist of nearly straight rolls. If the patches are targets or spirals, ^, as we have 
defined it, will be systematically smaller than i [Q. This is a consequence of basing the 
definition of ^ on the choice of the Fourier transform {i.e. on the standard definition of the 
structure factor), a transform based on a kernel consisting of plane waves. If the subpatterns 
have broadband structure in their Fourier transforms, this contributes to a broadening of the 
overall structure factor and a reduction of ^. The choice of the Fourier transform, which is 
sanctioned by long usage in condensed matter physics, is a natural one for the interpretation 
of scattering experiments, but does not correspond very well to what the eye naturally picks 
out as "coherent structures" in SDC patterns. Obviously, many other sorts of transform 
could be applied. Also, while the eye is drawn to the spirals, it is well to remember that 
SDC patterns contain significant areas which are dominated by other types of defects, such 
as grain boundaries, dislocations etc. The measured correlation length must be interpreted 
to be an average over of all of these, weighted according to their individual effect on the 
width of the structure factor. Similar considerations apply to the higher moments, to which 
we now turn. 

From the higher moments of S{k), we calculated the skewness s = /i3/i2^^^ and the 
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excess kurtosis k = /i4/i2 ^ — 3. The latter is defined in such a way that k is exactly zero 
for a gaussian distribution. Figure ^ shows the dependence of these statistics on e. For low 
e, the distribution is symmetric and indistinguishable from a gaussian, so that both s and 
K are nearly zero. The appearance of spirals near e ~ 0.25 is accompanied by a sudden 
increase in both statistics. After a considerable transient, the skewness continues to increase 
as the SDC develops. The excess kurtosis shows a very sharp jump at the onset of SDC, 
after which it is remarkably constant. These are nothing more than descriptive statistics, 
but they do contain the essential global features of the time-averaged spatial correlations of 
SDC that can be extracted in a model-independent way. The skewness has previously been 
used as a diagnostic of the onset of SDC 1^,0]. Our data indicates that the excess kurtosis 
gives an even clearer indication. Qualitatively, the increases at the threshold of SDC can be 
understood from the considerations of the previous paragraph. Spirals represent localized 
structures whose Fourier transforms contain significant asymmetric contributions in wings 
distributed about the local mean k. As the spirals appear and proliferate, S{k) develops 
both skewness and kurtosis. The larger scatter in the skewness for 0.2 < e < 0.5 relative to 
that of the kurtosis apparently reflects the fact that the fourth moment is less sensitive to 
sampling statistics in the calculation of S{k). In the lower part of this e range, the spirals 
appear only rather intermittently and the statistics of their contribution to the skewness is 
not very good. The extraction of moments is simply an objective way of disentangling the 
changes in the shape of S{k). The moments, particularly the higher ones, lack any simple, 
direct physical interpretation. However, such statistics may be useful for quantitatively 



comparing our results to numerical simulations |32|. 



D. Lifetimes from the three-dimensional structure factor 

We now return to the analysis of the structure factor S{k,uj), extracted from spacetime 
data sampled fast enough that successive pictures are correlated. Some insight into the de- 
gree of temporal correlation can be obtained by plotting two-dimensional slices of spacetime, 
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as in Fig. ^. The temporal behaviour of the pattern as it passes a particular line in space 
is shown spread out along the time axis. The regions of narrow stripes correspond to rolls 
which happen to meet the chosen line at nearly right angles, while the broad patches are 
rolls nearly parallel to the chosen line. A stack of 256 pictures like Fig. |^ make up the three 
dimensional spacetime data set at each e. Each data set represents 2^^ bits of information. 
The calculation of S{k, u), which destroys the phase information, reduces the number of bits 
by a factor of two. In addition, we ignored regions of k — u space far from the ring, which 
reduced the amount of data substantially. 

To reduce the data still further, we make use of the rotational symmetry of S{k, uj) about 
the k = axis, and take an azimuthal average to find S{k,uj). An example of the resulting 



two-dimensional data is shown in Fig. |T0[ In the k — uo plane, S decreases monotonically in 
the direction of increasing uo for all k. Because of the relatively rapid sampling, we cannot 
analyze as large a volume of spacetime as was the case for the data discussed in the previous 
section. As a consequence, statistical fluctuations are larger. 

The rate of decrease of S{k, uj) as a function of is a quantitative measure of the lifetime, 
or correlation time, of structures in the pattern with characteristic size l-n jk. Unfortunately, 
sampling statistics limit the useful part of the data to k values near (fc), i.e. near the peak 



of S{k,uj). Figure |TT] shows a plot of S{{k),uj) vs. for e = 0.75. If correlations decrease 



exponentially in time, this relation will be Lorentzian. The solid line shows a fit to 

^ = A/(l+cuV), (10) 

with the correlation time r and arbitrary amplitude A as fit parameters. In many cases, 
however, an exponential provides a better fit, and this also yields an arbitrary amplitude and 
a time scale as parameters. The times extracted from exponential fits were proportional to 
the correlation times found from Lorentzian fits. We adopted Lorentzian fits since they had 
the simplest interpretation in the time domain. In order to improve statistics, we averaged 
the lifetimes found for a range of k centered on (k); r is only a rather weak function of k. 
The dependence of r on e is shown on logarithmic scales in Fig. As with ^(e), we find 



a very rapid decrease of r with increasing e. A fit to a power law is shown by the straight 
fine. It gives an exponent of 1.43 ± 0.05. Thus the data are consistent with a threshold at 
e = and an exponent of about 3/2. 

The powerlaw analysis given above suggests a time scale which diverges at e = as 
1/e'^/^, whereas critical slowing down near bifurcation points usually involves a time scale 
proportional to 1/e. In order to see whether a 1/e dependence of r might also be consistent 
with the data, we show {r/ty)"^ vs. e in Fig. |13|. Also shown, as a dashed line, is the 
powerlaw fit with the threshold at e = 0. One sees that the data can be described equally 



well by a straight line in Fig. |13|, but this line passes through a shifted threshold at a 
positive Cc — 0.22. This value of ec is close to the onset of SDC given by other diagnostics 
as shown in Fig. |[ This interpretation of the data for r suggests that SDC appears via 
a bifurcation beyond which the abundance of spiral excitations grows algebraically. Other 
models which have been suggested imply an exponential growth of the frequency of spiral 
occurrence. At this time we know of no theoretical explanation of the powerlaw behaviour 
with a threshold either at e = or at Cc. It is worth noting that the data for ^ shown in 
Fig. 1^ can not be fit by a function with a non-zero threshold. 

IV. DISCUSSION 

The best established general conjecture about spatio-temporal chaos concerns the "ex- 
tensive" nature of the fractal dimension of the attractor which governs the spatio-temporal 
dynamics. [0,0] The dimension D is conjectured to be related to the characteristic size of 



the system L as 

D = {LIQ' (11) 

for "large" systems with L ^ C,c, where C,c is defined to be the "chaos" correlation length. 
Here d, the spatial dimension of the system, is 2 for our convection patterns. Indeed, the 
validity of Eq. |TI| might be considered as a defining feature of STC. Parts of a large system 
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are presumed to be nearly dynamically independent if they are separated by a distance of 
order C,c- The dimension D is extensive in the sense that it scales with the volume of the 
system, D oc L'^. Each volume ^c"^ acts like a dynamical subsystem (whose own dimension 
need not be small) which is weakly coupled to the whole, and contributes additively to the 
overall dimension D. This scaling scenario has been confirmed in a number of numerical 



simulations [p],^3 H4^ i |46[| in which D can be calculated as a function of L directly for certain 
model dynamical equations. One can then define an intensive dimension density 6 = D/L'^ 
so that 5 = 1/^/. 

The important question for the analysis of experiments is the relation, if any, between 
^c, the range over which the dynamics is physically correlated, and the observable pattern 
correlation length ^ as defined by Eq. p. If a simple proportionality (with a coefficient of 
order one) exists, we can estimate the rather high dimension of the attractor for SDC in a 
large but finite cell, a task which is essentially impossible by the usual time series methods 



used for low-dimensional systems. Encouragingly, several large scale simulations P2|- |4^J^ 
have found simple relations between C,c and ^. Simulations of the complex Ginzburg-Landau 
(CGL) equation in chaotic regimes P^-^ have found that the correlation length of the field 
amplitude is simply proportional to ^c- Within the parameter space of the CGL equation, 
however, one has transitions between phase- and defect-dominated chaotic states, and the 
overall phenomenology is not simple |4^j45[] . In general, one must choose carefully among 
the various possible definitions of C, for the proportionality of C,c and ^ to hold. For example, 



in simulations of the CGL equation [^^-Q , it was found that neither the correlation length 



derived from the complex field itself, nor that of its phase alone scaled linearly with C,c', only 
the correlation length derived from its amplitude did. A similar proportionality between 
and a suitably defined ^ has also been found in a completely different coupled-map lattice 
system Again, in this system other correlation lengths were found that did not scale 
linearly with C,c- There seems to be no prescription, other than a direct numerical test, for 
deciding which correlations lengths scale and which do not. 

Under the optimistic assumption that ^ ~ C,c, Eq- O with L = T and our result that 
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^ ~ e~^/^, imply that the dimension D for SDC increases hnearly with e, reaching 1400 at 
e = 1.4. That is, at e = 1.4, 1400 correlation areas vr,^^ fit in the area of the cell vrF^. The 
dimension density 6, which would also be linear in e, reaches ~ 0.25. This is comparable to 



that found in simulations of various model equations [^^cI'eII- course, this analysis is 
predicated on the untested assumption that ^ ~ ^c, which at present could only be estab- 
lished by prodigious amounts of numerical effort. In principle, the relation could be tested 



for SDC by finding the Lyapounov dimension D for realistic Navier-Stokes simulations [32 



or, less ambitiously, for a modified Swift-Hohenberg model [p9| -|31[|, using the direct meth- 
ods of Ref. Such a calculation, although very computationally expensive, is probably 



feasible [47 



The correlation lengths and times of SDC at large e must eventually be bounded by the 
cutoffs imposed by the basic scales of the convection process itself, the cell height d and 
the vertical time scale t^. At e ^ 1.4, we found C, ^ "2(1 and r ^ 6t„. Thus, the power 
law dependence on e cannot persist for much larger e. On the other hand, for small e, 
the experimental pattern made a crossover between SDC and the frustrated state, which is 
dominated by sidewall effects. A detailed study of the T dependence [§ of the onset of SDC 
would be most interesting. A very large F might be achieved by working close to the critical 
point of the gas [|^J28| J5|], but ultimately the experimental limits would be set by practical 
constraints on the acquisition and analysis of the very large amount of data that would be 
required. 



V. CONCLUSION 

In this paper we used a simple structure-factor analysis to characterize the space-time 
behaviour of shadowgraph images of SDC. The spatial and temporal correlations were ob- 
tained from the first few moments of an appropriately averaged structure factor. The length 
and time scales were found to scale as power laws in the reduced control parameter e. At 
present these power laws have no theoretical explanation. It is conjectured that the pattern 
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correlation-length is a measure of the size of physically correlated patches in the underlying 
chaotic dynamics. If this is so, the extensive dimension and the intensive dimension density 
of the dynamics both scale linearly with e. The data for the correlation time, although not 
definitive, suggest that there may be a sharp onset of SDC at €c > 0. Beyond these heuristic 
considerations, we have only descriptive statistics to characterize the general features of the 
motion. The higher moments of the structure factor which we have measured do not lend 
themselves to direct interpretation, but may be useful in the comparison of simulations to 
the data. In general, SDC serves to highhght the difficult challenge that extensive chaos 
presents to our understanding of dynamical systems. 
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FIGURES 

FIG. 1. The sequence of patterns observed as e is increased, (a): Nearly straight rolls at 
e = 0.043. (b): Weakly time dependent "frustrated" state at e = 0.245, just below the onset of 
SDC. 

FIG. 2. The sequence of patterns observed as e is increased, (a): At e = 0.536, above the onset 
of SDC. Spirals coexist with focus singularities at the edge of the cell, (b): At e = 0.894, fully 
developed SDC fills the cell. 

FIG. 3. The pattern observed at low e = 0.13, with the sidewall heater turned on. With the 
sidewall heater turned off, a pattern similar to Fig. P) is found instead. The heater has the effect 
of changing the roll orientation from perpendicular to parallel at the sidewall. 

FIG. 4. The patterns in the SDC regime at e = 0.644, (a) with the sidewall heater turned off, 
and (b) with the sidewall heater on. SDC is insensitive to the sidewall conditions, except very near 
the edges. 

FIG. 5. Patterns observed in the square cell with sidewall heating applied to the upper and 
lower edges of the ceh. (a) e = 0.17, showing a "frustrated" state, (b) SDC at e = 0.77. 
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FIG. 6. The time average of kS{k) for e = 0.043 (round points), e = 0.306 (triangles), and 
e = 1.032 (squares). 
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FIG. 7. A plot of 1/^^ vs. e. The solid line corresponds to ^ = 0.194e. The good fit to the 
data demonstrates that ^ ~ e~^^'^ for small e. 
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FIG. 8. The higher moments of the time averaged structure factor, (a) the skewness and (b) 
the excess kurtosis, as a function of e. 
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FIG. 9. A snapshot of the spacetime structure of fully developed SDC in the circular cell, (a) 
A space-space shadowgraph image at e = 1.032. The image is 106d wide, (b) A spacetime slice 
of SDC, with the space coordinate horizontal and time running downward. The image represents 
spacetime data along a horizontal line near the midpoint of image (a), spread out over a time 
interval of 363t^ . The dimensions of (b) correspond to 46^ in the space direction and 55r in time. 
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FIG. 10. A plot of the surface of a typical azimuthally averaged 5(11^1, a;), showing the relation 
between the various statistical quantities. 
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FIG. 11. A dimensionless plot of S{k,uj) vs Lot^, for k = (k), at e = 0.747. The solid line is a 
fit to a two parameter Lorentzian. The fit gives r = 12.3t„. 
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FIG. 12. The correlation time r extracted from Lorentzian fits like the one shown in Fig. 11 



vs. e, for /c in a narrow range around (k). The soHd line is a fit which gives an exponent of 
-1.43 ±0.05. 
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FIG. 13. A plot of vs. e, for the same data as in Fig. 12. This shows that r ~ (e — ec)~^ 
with Cc — 0.22 is consistent with the data. The powerlaw fit with an exponent of 3/2 and a 
threshold at e = is given by the dashed line. 
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